Is the Kelvin Theorem Valid for High-Reynolds-Number Turbulence? 
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The Kelvin-Helmholtz theorem on conservation of circulations is supposed to hold for ideal in- 
viscid fluids and is believed to be play a crucial role in turbulent phenomena, such as production of 
dissipation by vortex line-stretching. However, this expectation does not take into account singular- 
ities in turbulent velocity fields at infinite Reynolds number. We present evidence from numerical 
simulations for the breakdown of the classical Kelvin theorem in the three-dimensional turbulent 
energy cascade. Although violated in individual realizations, we find that circulations are still con- 
served in some average sense. For comparison, we show that Kelvin's theorem holds for individual 
realizations in the two-dimensional enstrophy cascade, in agreement with theory. The turbulent 
"cascade of circulations" is shown to be a classical analogue of phase-slip due to quantized vortices 
in superfluids and various applications in geophysics and astrophysics are outlined. 
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The theorem on conservation of circulations due to 
Helmholtz [1] and Kelvin [2] is a fundamental fluid dy- 
namical result with many important consequences. The 
circulation invariants are a topological obstruction to 
ideal vortex-reconnection in classical and quantum flu- 
ids [3,4]. Ertel's theorem [5] on conservation of poten- 
tial vorticity in geophysical fluid dynamics is a differen- 
tial version of the Kelvin theorem [6]. The conservation 
of circulations was also argued by G. I. Taylor [7-9] to 
play a key role in the enhanced production of dissipa- 
tion in turbulent fluids, by the process of vortex line- 
stretching. More recently, regularizations of the Navier- 
Stokes equation have been proposed as model equations 
for large-scale turbulence [10,11], motivated by requir- 
ing that a Kelvin theorem be preserved. However, de- 
spite the frequent application to turbulent flows, existing 
proofs of conservation of circulation are valid only for 
smooth, laminar solutions of the ideal fluid equations. 
One might naively expect that the conservation will hold 
better as the viscosity decreases, or the Reynolds num- 
ber increases, but this has never been shown. This Let- 
ter presents evidence from numerical simulations both for 
the breakdown of the classical Kelvin theorem and for a 
generalized version that is still valid in turbulent flow at 
infinite Reynolds number. 

We consider the incompressible Navier-Stokes equation 
in space dimension d > 2 : 



d t u + (u-V)u = -Vp + v A u, V-u = 0, 



(1) 



where u(x, t) is the velocity, p(x, f) the kinematic pres- 
sure (or enthalpy), and v the kinematic viscosity. The 
classical Kelvin theorem in this context states that for 
any closed, rectifiable loop C at an initial time to, the 
circulation T{C,t) — § c u\ u(t)-dx satisfies 



Au(i)-dx, 



(2) 



where C(t) is the loop advected by the fluid velocity, at 
time t. In the inviscid limit — > 0, it follows heuristically 
from equation (2) that the circulation is conserved for 
any initial loop C. This conservation law can be shown 
to be a consequence of Noether's theorem for an infi- 
nite dimensional gauge symmetry group of the ideal fluid 
model associated to particle-relabelling [12,13]. 

The conservation can be anomalous, however, in the 
zero- viscosity limit, if the velocity field becomes too sin- 
gular. In that case, the righthand side of (2) need not 
vanish as v — > 0. The situation is quite similar to that 
for conservation of energy, which, as observed by On- 
sager [14], can also become anomalous in the sense that 
the energy dissipation e = v\ Vu| 2 need not vanish in the 
inviscid limit. For the energy integral, the existence of a 
dissipative anomaly can be associated to a turbulent cas- 
cade with a constant mean flux of energy to arbitrarily 
small length-scales. See [15-17]. 

Anomalous circulation conservation may be formu- 
lated in a similar fashion [18]. Let = Gi * u de- 
note the low-pass filtered velocity at length-scale £, where 
Ge (r) = £~ d G(r/£) is a smooth filter kernel. Then u ( sat- 
isfies the following equation (dropping viscous terms that 
are small in the inertial-range) : 



d t \it + (u r V)u e = -Vp t + i e 



(3) 



where p e is the large-scale (modified) pressure and where 
is the turbulent vortex-force: 



te = (uxu), — UfXU( 



(4) 



Let us choose a rectifiable closed loop C in space. We 
define Ci(t) as the loop C advected by the filtered veloc- 
ity Tit and a "large-scale circulation" with initial loop C 
as the line-integral Ti(C,t) = ^ , t -.viz(t) ■ dx. It follows 

from (3) that 



(d/dt)T t {C, t) 



C(t) 



U(t) ■ dx 



(5) 



C e (t) 



1 



2 



The righthand side of (5) represents a turbulent trans- 
port of lines of large-scale vorticity ZJe(t) out of the loop 
Ci(t). Using formula (4) and the righthand rule, it is 
easy to see that fe provides a torque that reduces the 
large-scale circulation around the loop when a vortex-line 
migrates outside it. This process can be seen to be local- 
in-scale, with most of the transport due to scales close 
to £. Therefore, it is natural to think of the process as a 
"circulation-cascade" [18]. This motivates the definition, 
for any loop C and filter length £, of a circulation-flux 
Ke(C) = — § c fe ■ dx. If the effect of the modes at lengths 
< £ is, on average, to diffuse the lines of large-scale vortic- 
ity, then this definition implies that the signs of Tg(C, t) 
and Ke(C, t) should be positively correlated. 

The righthand side of (5) — or the circulation-flux — 
does not need to vanish in the limit taking first v — > 0, 
then £ — > 0. This provides a purely inviscid mechanism 
for violation of the Kelvin theorem, if the fluid velocity 
is sufficiently singular. A simple estimate of the vor- 
tex force shows that |f)| = 0(\Su(£)\ 2 /£), where 5u(£) is 
the velocity increment across a length £ [18]. Suppose 
that the velocity is Holder continuous in space with ex- 
ponent h, so that |(5u(€) | = 0(£ h ). It follows then that 
K t {C) = O (£ 2h - 1 ) for any loop C of finite length [18]. 
Thus, the Kelvin theorem will be valid under these as- 
sumptions if h > 1/2. The latter condition will hold in 
certain cases, e.g. the 2c? cnstrophy cascade, where it is 
expected that h — 1 (with logarithmic corrections). How- 
ever, in the 3d energy cascade, the condition is not ex- 
pected to be true, since even the mean-field Kolmogorov 
exponent ft = l/3is<l/2. Thus, the rigorous results do 
not settle the issue of whether Kelvin's theorem is valid 
for infinite Reynolds-number turbulence in 3d, where a 
non-vanishing circulation-flux is possible. 

To explore these questions, we have carried out a direct 
numerical simulation (DNS) of eq. (1) on a 1024 3 peri- 
odic grid using a pseudo-spectral parallel code with full 
dealiasing and time-stepping by a second-order Adam- 
Bashforth method. The kinetic energy was forced in 
the first two shells [19] and a statistical stationary state 
was achieved after 5 large-eddy turnover times. The fi- 
nal Taylor-scale Reynolds number was Re x — 383 and 
about a decade of incrtial range was achieved with con- 
stant mean energy flux and a spectral exponent close to 
—5/3. We investigated the statistics of the circulation 
and circulation-flux on square loops of edge-length R. 
The line-integrals were evaluated by using Stokes theo- 
rem and then calculating the surface-integrals over the 
square faces. For each value of R all square loops were 
considered on the computational grid with that edge- 
length, at all positions in the periodic domain and with 
all three possible orientations. 

In Fig. 1(a) plot the probability density function 
(PDF) of the circulation-flux Kg(C) for a square loop 
with radius R = 64 (lattice units) and for several filter 
lengths I < R. As £ decreases through the inertial-rangc, 
the PDF's approach a form that is nearly independent 
of the filter length. This invariant PDF is non-Gaussian 



with stretched exponential tails, qualitatively similar to 
the PDF of the circulation itself [20]. The crucial point 
here is lack of dependence of the PDF on the filter-length 
£, which is analogous to the independence of the mean en- 
ergy flux in the inertial-range. Note that the mean of the 
circulation-flux must be zero in homogeneous turbulence 
because of the exact identity fg = — V-t^ + Vk^ for the 
turbulent vortex-force, where Tg = (uu)< — u^u^ is the 
turbulent subscale stress and Kg — (l/2)trr£ is the sub- 
scale kinetic energy. However, higher-order moments of 
Kt(C) can be non- vanishing. In Fig. 1(b) we plot the rms 
value of the circulation-flux as a function of k c = w/£ for 
several values of R in the inertial range. The near inde- 
pendence of rms flux to k c verifies the existence of the cir- 
culation cascade. In an infinitely long inertial-range, for 
v — > 0, this nonlinear cascade would lead to breakdown 
of the classical Kelvin theorem. It is also interesting to 
consider the i?-dependence of the inertial-range value of 
the rms circulation- flux. In the inset to Fig. 1(b) we plot 
the rms value versus R, revealing a scaling very close to 
i? 1 / 2 . Since Kg(C) has the dimension of velocity-squared, 
the mean-field Kolmogorov scaling would be i? 2 / 3 . The 
smaller exponent shows that rms circulation-flux scales 
anomalously with R, similar to the circulation itself [20] . 




FIG. 1. (a) PDF of the circulation-flux for loops with ra- 
dius R = 64 and for several filter lengths £ < R. (b) The rms 
value of the circulation-flux as a function of k c for various 
loop sizes R. The inset plots the plateau rms value versus R. 
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It is illuminating to compare these numerical results for 
3d with corresponding results for the enstrophy cascade 
in 2d. We have analyzed the solutions of a 2048 2 DNS, 
the details of which are given in [21]. This simulation 
yielded about a decade and a half of inertial-range with 
constant mean enstrophy-flux and an energy spectral ex- 
ponent close to —3. In Fig. 2(a) we plot the PDF's of 
circulation-flux Ke(C) from this DNS for a square loop 
with R = 128 and several filter-lengths £ < R. Unlike 
the corresponding results in Fig. 1(a) for 3d, the PDF's 
are not independent of £ but instead narrow rapidly as 
£ decreases. This is quantified in Fig. 2(b), which plots 
the rms circulation- flux versus k c , again for loops with 
radius R = 128. A power-law decay is observed for in- 
creasing k c with an exponent between —1 and —2, con- 
sistent with the rigorous bound in [18]. These results 
show that the Kelvin theorem is valid in the 2d enstro- 
phy cascade, whereas circulations are not conserved in 
the 3d energy cascade range as £ — > 0, due to persistent 
nonlinear transport of vortex-lines. 




FIG. 2. (a) The PDF's of circulation-flux from a 2048 2 sim- 
ulation of 2D enstrophy cascade for loops of radius R = 128 
and several k c = ir/£. (b) The rms circulation-flux plotted 
versus k c , also for R = 128. 

Although the Kelvin theorem is not valid in the clas- 
sical sense in 3d turbulence, it was conjectured in [18] 
that a weaker form may still be valid. Consistent with 
the results above, the circulation T(C, t) on an advected 
loop C(t) is expected not to be invariant in time. In 



fact, associated with the phenomenon of "spontaneous 
stochasticity" in the zero- viscosity limit [22,23], the loop 
C (t) should be random for a fixed initial loop C and ad- 
verting velocity field u. Under these circumstances, the 
time series of F(C, t) will be a stochastic process. It was 
proposed in [18] that the random time-series of circula- 
tions should possess the "martingale property" : 

(T(C,t)\T(C,s), s<t')=T{C,t'), t>t'. (6) 

That is, the conditional expectation of the circulation 
in the future should be the last given value. This is a 
natural statistical generalization of the Kelvin theorem 
and, in fact, is a generalized form of the inviscid Eulcr 
equations of motion for a turbulent fluid [18]. 

A full test of these ideas will be quite difficult and in- 
volve, among other things, careful Lagrangian tracking 
of the loops C(t). Here we check a somewhat weaker re- 
sult. In Fig. 3 we plot the conditional average circulation 
flux (K e (C)\T(C) = T) from the 3d DNS as a function of 
k c = ir /I, for a square loop with R — 64 and for various 
values of the circulation level T. Unlike (6), this expec- 
tation includes an average over the turbulent ensemble 
of velocities u. The plot in Fig. 3 shows that the sign of 
Ki(C) is positively correlated with that of r(C), consis- 
tent with the diffusive character of turbulent vortex-line 
transport. Thus, the effect of subscale modes at length- 
scales < £ will be to reduce the magnitude of the circu- 
lation, regardless of its sign. However, the conditional 
average flux for each value of T tends to zero as £ de- 
creases through the inertial range. This is a true effect 
of the nonlinear dynamics, as illustrated by the dashed 
line in Fig. 3, which plots the expectation of the viscos- 
ity term v § c Aufdx conditioned upon T = T rms , as a 
function of k c . Clearly viscous effects are negligible over 
the range considered. These results show that the Kelvin 
theorem, although violated in individual realizations, is 
still valid for 3d turbulence in some average sense. 
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FIG. 3. The conditional average circulation-flux from the 
3d DNS as a function of k c , for a square loop with R = 64 
and for various values of the circulation level (solid lines) . The 
dashed line shows the conditional expectation of the viscosity 
term, given F = IY ms , as a function of k c . 
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It is worth pointing out that the "circulation cascade" 
verified for fluid turbulence in this work is a classical ana- 
logue of the "phase slip" phenomenon in superfluids, such 
as 4 He below the A-point [24]. For example, consider the 
decay of a superfluid flow in a thin toroidal ring [25] . The 
decay is mediated by the (thermal or quantum) nucle- 
ation of quantized vortices which migrate out of the ring. 
The passage of a vortex across the toroidal cross-section 
induces by phase-slip a pulse of torque which decreases 
the circulation around the ring. It was emphasized al- 
ready by Anderson [24] that phase-slip occurs in superflu- 
ids only because the quantized vortices are not material 
objects moving with the fluid. This is possible because 
of the singular vortex core, where the superfluid density 
drops to zero. The mechanism of turbulent circulation- 
cascade that we have considered is quite similar. The 
vortex lines in a turbulent flow are also not material, 
because singularities in the velocity field allow them to 
diffuse relative to the fluid. Unlike in superfluids, this 
is a continuous process, since classical vortices are not 
quantized. There is also no need for the singularities to 
be nucleated as fluctuations, since they are everywhere 
present in the turbulent flow. Finally, we note that in 
turbulence the diffusion of vortex-lines is not persistent 
in scale, on average, and does not lead to irreversible de- 
cay of circulations. As shown by Fig. 3, the partial effect 
of the modes adjacent in scale (£ ~ R) leads to a mean 
decay but, for £ -C R, the circulation is conserved on av- 
erage. In this limit, time-reversal symmetry is restored. 

The results in this work have many important impli- 
cations. The turbulent transport of vortex-lines, per- 
sistent as I — > 0, provides an inviscid mechanism for 
vortex-reconnection and other changes of line topology. 
However, a statistical form of the Kelvin theorem seems 
to survive, which is crucial to justify Taylor's vortex- 
stretching mechanism of turbulent energy dissipation 
[7-9]. Further research must clarify exactly to what ex- 
tent conservation of circulation survives for turbulent so- 
lutions of fluid equations. Numerical studies must be 
extended to Lagrangian tracking of advected loops C(t). 
Theoretical investigations are required to take into ac- 
count the presumably fractal nature of such loops [26]. 
Circulation-cascade should occur also in superfluid tur- 
bulence [4], mediated in part by phase-slip of quan- 
tized vortex-lines. Results similar to the present ones 
may be developed also for conservation of magnetic flux 
(Alfvcn's theorem [27]), important in magnetohydrody- 
namic turbulence, and Ertel's theorem [5] in geophysical 
turbulence. The breakdown of Alfven's theorem and the 
"frozen-in" character of magnetic flux-lines is of partic- 
ular importance to account for the fast reconnection ob- 
served in a variety of astrophysical turbulent flows [28] . 
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